TTP 95-06 

hep-ph/ 9503272 



INCLUSIVE DECAYS 
OF HEAVY FLAVOURS 



Thomas Mannel 

Institut fur Theoretische Teilchenphysik, University of Karlsruhe 
Kaiserstr. 12, D - 76128 Karlsruhe, Germany. 



Abstract 

Recent progress in the theoretical description of inclusive heavy flavour 
decays is reviewed. After an outline of the theoretical methods appli- 
cations to total decay rates and semileptonic decay spectra are pre- 
sented. 



Contribution to the workshop on "Heavy Quark Physics", 
December 13 - 16, 1994, Bad Honnef, Germany. 



TTP 95-06 

March 1995 



1 Introduction 



Decays of heavy flavours play an important role in the determination of the 
not yet well explored CKM sektor of the standard model. The weak pro- 
cesses that need to be studied involve transitions among quarks of different 
flavours; however, to connect the observed transitions among hadrons to the 
underlying quark processes one has to deal with the bound state problem of 
strong interactions. 

For the decays of heavy flavours one may take advantage of the fact that 
the mass rriQ of the heavy quark is large compared to the scale A, which is 
determined by the light degrees of freedom and thus is of the order of Aq C d- 
The amplitudes or the transition rates for the decays of heavy flavoured 
hadrons are evaluated as an expansion in powers of A/tuq, the leading term 
of which corresponds to an infinitely heavy, static quark p], @j ■ 

The 1 / rriQ expansion has been formulated in the laguage of effective field 
theory, the so called Heavy Quark Effective Theory (HQET) 0, [|. In this 
framework one may systematically access the mass dependence of any matrix 
element involving heavy hadron states and heavy quark fields. Furthermore, 
in the infinite mass limit two additional symmetries appear 0, which play 
an important role in the context of exclusive heavy flavour decays. 

In the past few years methods have been developed to apply the I/tuq 
expansion also to inclusive decays of heavy hadrons [|, |, ||, [7|, || |9|, [H| , and 
the purpose of this talk is to give a short review of these developments. It is 
impossible to cover all topics in the field of inclusive decays due to length and 
time restrictions and consequently this review can cover only some selected 
topics. 

For inclusive decays a 1/mq expansion is obtained for the rates by an 
approach similar to the one known from deep inelastic scattering. The first 
step consists of an Operator Product Expansion (OPE) which yields an in- 
finite sum of operators with increasing dimension. The dimensions of the 
operators are compensated by inverse powers of a large scale, which is in 
general of the order of the heavy mass scale. The decay probability is then 
given as forward matrix elements of these operators between the state of the 
decaying heavy hadron; these matrix elements still have a mass dependence, 
which then may be extracted in terms of a 1/rriQ expansion using HQET as 
for exclusive decays. 

In section 2 we give a short description of the OPE approach to inclusive 
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decays. Section 3 summarizes the results for the total rates, including the 
first non-trivial non-perturbative corrections. In section 3 the method is 
applied to the charged lepton energy spectrum in inclusive semileptonic B 
decays. It turns out that the endpoint region cannot be described by a 
1/rng expansion; rather a partial resummation of the I/tuq expansion is 
required, which is closely analogous to the leading twist term in deep inelastic 
scattering [11, O, 13, Il4j]. This is summarized in section 4. Finally, we 



conclude and point out a few open questions which are currently under study. 



2 Operator Product Expansion 

The effective Hamiltonian for a decay of a heavy (down-type) quark is in 
general linear in the decaying heavy flavoured quark 

We// = QR (i) 

where the operator R describes the decay products. In the following we shall 
consider semileptonic decays, for which 

R sl = ^=V Qq 7/l (l - l,)q (^7 M (1 - 75^), (2) 

where q is an up-type quark (c or u, since we shall consider b decays). Sim- 
ilarly, for nonleptonic decays the Cabbibo allowed contribution corresponds 
to 

Rm = ^F Q? y; ? ,4(C + (m 6 )+C'_(m 6 ))7 M (l-75)g(g / 7 M (l-75)? ,/ ) (3) 

+(C + (m b ) - C7_(m 6 )) 7M (l - IsW - 7 5 )<z)] , 

where q' (q") is a down-type (up-type) quark and Vq q the corresponding CKM 
matrix element. The coefficients C±(mb) are the QCD corrections obtained 
from the renormalization group running between M\y and m&; in leading 



logarithmic approximation these coefficients are |T5 

r2 

,2 



a s (M^) 6 1 

C±{mb) = [ ^KT ] ' Wlth 7+ = 3^2^ = "2 7 "' (4) 

where a s (fi) is the onle-loop expression for the running coupling coupling 
constant of QCD. 
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Finally, for radiative rare decays we have 

R rare = -^V tb V; s C 7 {m b )^m b a^(l + j 5 )sF^, (5) 

where C-j{m b ) is again a coefficient obtained from running between My/ and 
m b . Its value is C-?(rrib) ~ 0.3, the corresponding analytical expression may 
be found in ]16| . 

The inclusive decay rate for a heavy hadron H containing the quark Q 
may be related to a forward matrix element by 

T oc ^2(2-k)*5\P h - P x )\(X\H eff \H(v)}\ 2 (6) 
x 

= Jd 4 x(H(v)\H eff (x)Hl ff (0)\H(v)) 

= 2 Im / d 4 x (H(v)\T{H eff (x)Hi ff (0)}\H(v)), 



where \X) is the final state which is summed over to obtain the inclusive 
rate. 

The matrix element appearing in ([]) contains a large scale, namely the 
mass of the heavy quark. The first step towards a 1/rriQ expansion is to 
make this large scale explicit. This may be done by a phase redefinition 
which leads to 

r oc 2 Im/ d 4 xe- m ^ x (H(v)\T{n eff (x)Hl ff (0)}\H(v)), (7) 

where 

Heff = QvR Qv = e~ imQVX Q. (8) 

This relation exhibits the similarity between the cross-section calculation in 
deep inelastic scattering and the present approach to total rates. In deep in- 
elastic scattering there appears a large scale which is the momentum transfer 
to the leptons, while here the mass of the heavy quark appears as a large 
scale. 

The next step is to perform an operator product expansion of the product 
of the two Hamitonians. After the phase redefinition the remaining matrix 
element does not involve large momenta of the order of the heavy quark mass 
any more and hence a short-distance expansion becomes useful, if the mass 
rriQ is large compared to the scale A determining the matrix element. The 
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next step is thus to perform an operator-product expansion, which has the 
general form 



J d ^e^(H(v)\T{n eff (x)Hl ff (0)}\H(v)) (9) 

OO / -1 \ n 

= £ (^J (H(v)\O n+ 3\H(v))„ 

where O n are operators of dimension n, with their matrix elements renor- 
malized at scale /i, and C n are the corresponding Wilson coefficients. These 
coefficients encode the short distance physics related to the heavy quark 
mass scale and may be calculated in perturbation theory. All long distance 
contributions connected to the hadronic scale A are contained in the matrix 
elements of the operators O n+3 . 

Still the matrix elements of O n+ z are not independent of the heavy quark 
mass scale, but this mass dependence may be expanded in powers of 1/mg 
by means of heavy quark effective theory. This is achieved by expanding 
the heavy quark fields appearing in the operators O n as well as the states by 
including the corrections to the Lagrangian as time-ordered products. In this 
way the mass dependence of the total decay rate may be accessed completely 
within an expansion in 1/rriQ. 

The lowest-order term of the operator product expansion are the dimension- 
3 operators. Due to Lorentz invariance and parity there are only two com- 
binations which may appear, namely Q V ^Q V or Q V Q V . Note that the Q v 
operators differ from the full QCD operators only by a phase redefinition, 
and hence Q v j)Q v = Qj)Q and Q V Q V = QQ. The first combination is propor- 
tional to the Q-number current Q^^Q, which is normalized even in full QCD, 
while the second differs from the first one only by terms of order I/rug 
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QvQv = v^Qv^Qy + 77^2 K (iDf - (ivDf + % -a^G^ u h v + 0(l/m 6 Q ), 

Zniq L I . 

(10) 

where G^ u is the gluon field strength. 

Thus the matrix elements of the dimension-3 contribution is known to be 
normalized; in the standard normalization of the states this implies 

(H(v)\0 3 \H(v)) = (H(v)\Q v jQ v \H(v)) = 2m H , (11) 
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where m# is the mass of the heavy hadron. To lowest order in the heavy mass 
expansion we may furthermore replace tuh = ttlq and hence we may evaluate 
the leading term in the 1/mg expansion without any hadronic uncertainty. 
Generically the dimension-3 contribution yields the free quark decay rate. 
This has been previously used as a model for inclusive decays, but now it 
turns out to be the first term in a systematic 1/rrig expansion of total rates. 

A dimension-4 operators contains an additional covariant derivative, and 
thus one has matrix elements of the type 

(H(v)\0 4 \H(v)} cc {H(v)\Q v TD^Q v \H(v)) = A T v^. (12) 

Since the equations of motion apply for this tree level matrix element, one 
finds that the constant At has to vanish, and thus there are no dimension- 
4 contributions. This statement is completely equivalent to Lukes theorem 
[ Q| , since we are considering a forward matrix element, i.e. a matrix element 
at zero recoil 



The first non-trivial non-perturbative contribution comes from dimension- 
5 operators and are of order 1/toq. For mesonic decays there are only the 
two parameters Ai and A2 corresponding to matrix elements involving higher 
order terms that appear in the effective theory Lagrangian 

_ (H{v)\h v {iDyh v \H{v)) 

Al " 2M H (13) 

(H(v)\h v a^iDWh v \H(v)) 

Aa " 2M H ' (M) 

where the normalization of the states is chosen to be (H(v)\h v h v \H(v)) = 
2Mh, where Mjj is the mass of the heavy meson in the static limit. These 
parameters may be interpreted as the expectation value of the kinetic energy 
of the heavy quark and its energy due to the chromomagnetic moment of the 
heavy quark inside the heavy meson respectively. 

The parameter A2 is easy to access, since it is related to the mass splitting 
between H(v) and H*(v,e). From the 5-meson system we obtain 

A 2 (m 6 ) = ^(Mff* - M H ) = 0.12 GeV 2 ; (15) 

from the charm system the same value is obtained. This shows that indeed 
the spin-symmetry partners are degenerate in the infinite mass limit and the 
splitting between them scales as 1/mg. 
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The parameter Ai appearing is not simply related to the hadron spectrum; 
from the definition of Ai one is led to assume Ai < 0; a more restrictive 
inequality 

-Ai>3A 2 (16) 

has been derived in a quantum mechanical framework in |TJ] and using heavy- 
flavour sum rules ||19|| . Furthermore, there exists also a QCD sum rule esti- 
mate [20 1 for this parameter: 

Ai = -0.52 ±0.12 GeV 2 . (17) 



3 Total Decay Rates 



In this section we collect the results for the total rates including the first 
non-trivial non-perturbative correction. 

Inserting R s i as given in (0) one obtains for the total inclusive semileptonic 
decay rate B — > X c tv 



r 



B^X c lv 



G\m\ 

1927T 3 



cb\ 



1 + 



XT 
2m? 



9A 2 
2ml ' 



rric 



where the two /,■ are phase-space functions 



h{x) 



1 - 8s 2 + 8s 6 - a; 8 - 24x 4 logx, 

8 5 

1 x 2 - 8x 4 + 8x 6 + -a; 8 + 8a; 4 log x. 

3 3 



The result for B — > X u £u£ is obtained from (|T8| ) as the limit m c 
replacement V c b — > V ub 



(18) 



(19) 



and the 



B^X u lv 



G 2 F ml 

1927T 3 



ub | 



1 + 



Ai - 9A 2 
2ml 



(20) 



As was discussed above, the leading non-perturbative corrections in (|T8|) and 
( |2T)| ) are parametrized by Ai and A 2 . Estimates for these parameters have 
been discussed in section 2; in order to estimate the total effect of the non- 
perturbative effects we insert a range of values —0.3 > Ai > —0.6 GeV 2 ; 
from this we obtain 

^~-<S...4)* (21) 
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This means that the non-perturbative contributions are small, in particular 
compared to the perturbative ones, which have been calculated some time 
ago [^j], |22J . For the decay B — > X u lvt the lowest order QCD corrections are 
given by 



1927T 3 



w, 



ub 



1 + 



2a /25 



3vr V 4 



7T 



0.851V, 



ub 



1 h 



(22) 



and thus the typical size of QCD radiative corrections is of the order of ten 
to twenty percent. 

Similarly, one obtains the result for non-leptonic decays as 



' 1927T 3 



\m b J 



1 + t-^(Ai -9A 2 



2m\ 



(23) 

where the coefficients A, are expressed as combination of the Wilson coeffi- 
cents C±{m b ) given in (|) 



A 1 = -[C 2 (m b )+2Cl(m b )}, 



~[Cl(m b ) - C 2 _(m b )}, (24) 



and f 3 (x) 



1 



x 



2\3 



is another phase space function. Again the non- 



perturbative corrections turn out to be small, in the region of a few percent 
compared to the leading term, and the perturbative corrections turn out to 
be much larger than this. 

Finally, for the rare decay B — > X s ^ one may as well calculate the non- 
perturbative contribution in terms of Ai and A2. One obtains 



aG 



16tt 4 



W b \V ts V: d \ 2 \C 7 {m b )f 



1 + 



1 



2m\ 



(Ai-9A 2 ) 



(25) 



and the relative size of the nonperturbative corrections is the same as in the 
B — > X u lvi decays. 

Typically the non-pertubative corrections are much smaller than the ra- 
diative corrections. The only exception is the endpoint region of lepton 
energy spectra which receives both large perturbative as well as nonpertur- 
bative corrections. However, this is only a small region in phase space and 
the corrections to the total rates remain moderate. 
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4 Lepton Energy Spectra 



The method of the operator-product expansion may also be used to obtain 
the non-perturbative corrections to the charged lepton energy spectrum. In 
this case the operator product expansion is applied not to the full effective 
Hamiltonian, but rather only to the hadronic currents. The rate is written 
as a product of the hadronic and leptonic tensor 



dr 



4m 



B 



\V Qq \ 2 W^d(PS), 



(26) 



where d(PS) is the phase-space differential. The short- distance expansion is 
then performed for the two currents appearing in the hadronic tensor. Re- 
defining the phase of the heavy-quark fields as in (§) one finds that the mo- 
mentum transfer variable relevant for the short-distance expansion is mgv—q, 
where q is the momentum transfer to the leptons. 

The structure of the expansion for the spectrum is identical to the one 
of the total rate. The contribution of the dimension-3 operators yields the 
free-quark decay spectrum, there are no contributions from dimension-4 op- 
erators, and the I /ml corrections are parametrized in terms of Ai and A 2 . 
Calculating the spectrum for B — > X c tv yields || 0, |3|. 



dY 

dy 



G 2 F \V cb \ 2 


m\ 


1927T 3 




+ A1 




[m b (l - 




3A 2 




mf(l - 





9(1 - y - p)y 2 [{3(1 - p){\ - R 2 ) - 2y{\ - R 3 )} (27) 
Ai 



: (3R 2 -4R S 



ml(l-y) 



{R 2 - 2R 6 



3ml 



A 



+— 1[(6 + by) - 12R - (9 - 5p)R 2 + 10R d 



mi 



+ (A/[m b (l-y)}Y 



where we have defined 



and 



,m b J l-y 



y = 2E e / m b 



(28) 
(29) 
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Figure 1: The electron spectrum for free quark b —>■ c decay (dashed line), 
free quark b — > u decay (grey line), and B — > X c ev e decay including I /ml 
corrections (solid line) with Ai = —0.5 GeV 2 and A2 = 0.12 GeV 2 . The figure 
is from ||. 



is the rescaled energy of the charged lepton. 

This expression is somewhat complicated, but it simplifies for the de- 
cay B — > X u iv since then the mass of the quark in the final state may be 
neglected. One finds 



dT 

dy 



V ub \ 2 ml 



192vr 3 
Ai + 33A 2 



(2, 2 (3-2,) + ^4 + 2,(6 + 5,)4 N 



6(1 



3m 2 b 



5(1-2/) 



Ai 

3m 2 b 



-y) 

(30) 



Figure |I] shows the distributions for inclusive semileptonic decays of B 
mesons. The spectrum close to the endpoint, where the lepton energy be- 
comes maximal, exhibits a sharp spike as y — > y m ax- In this region we have 



_«9(l-y-p) 



2 + 



Ai 



P 



(m Q (l-y)) 2 \1 -p 




(31) 



which behaves like 5-functions and its derivatives as p — > 0, which can be seen 
in (pOf). This behaviour indicates a breakdown of the operator product ex- 
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pansion close to the endpoint, since for the spectra the expansion parameter 
is not 1/m.Q, but rather l/(mg — qv), which becomes l/(m<g[l — y}) after the 
integration over the neutrino momentum. In order to obtain a description of 
the endpoint region, one has to perform some resummation of the operator 
product expansion. 



5 Resummation in the Endpoint Region 

Very close to the endpoint of the inclusive semileptonic decay spectra only 
a few resonances contribute. In this resonance region one cannot expect to 
have a good description of the spectrum using an approach based on parton- 
hadron duality; here a sum over a few resonances will be appropriate. 

In the variable y the size of this resonance region is, however, of the order 
of (A/ttiq) 2 and thus small. In a larger region of the order A/itlq, which we 
shall call the endpoint region, many resonances contribute and one may hope 
to describe the spectrum in this region using parton-hadron duality. 

It has been argued in |Tl| that the 5-function-like singularities appearing 



in (|30D may be reinterpreted as the expansion of a non-perturbative function 
describing the spectrum in the endpoint region. Keeping only the singular 
terms of (^(J) we write 

^ = 2y\3-2y)S(y), (32) 
T b dy 

where 

oo 

S(y) = &(l-y) + ^a n 5^(l-y) (33) 

n=0 

is a non-perturbative function given in terms of the moments a n of the spec- 
trum, taken over the endpoint region. These moments themselves have an 
expansion in 1/mq such that a n ~ l/mg +1 , and we shall consider only the 
leading term in the expansion of the moments, corresponding to the most 
singular contribution to the endpoint region. 

Comparing (|30"D with (|32|) and ( |33"D one obtains that 

a = Jdy(S(y)-&(l-y)) = 0(l/ml) } (34) 
ai = Jy(S(y)-e(l-y)) = -^ F , (35) 
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where the integral extends over the endpoint region. 

The non-perturbative function implements a resummation of the most 
singular terms contributing to the endpoint and, in the language of deep 
inelastic scattering, corresponds to the leading twist contribution. A similar 
approach based on the parton model will be described in a separate talk at 
this conference 
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This resummation has been studied in QCD HT2| , [13| and the function 
S(y) may be related to the distribution of the light cone component of the 
heavy quark residual momentum inside the heavy meson. The latter is a 
fundamental function for inclusive heavy-to-light transitions, which has been 
defined in 







i 



2m 



B 



-{H{v)\h v 6{k+-iD + )h v \H{v)), 



(36) 



where k + = k + fc 3 is the positive light cone component of the residual 
momentum k. The relation between the two functions S and / is given by 



A 

S(y) = — I dk + f(k + ), 



(37) 



-m,Q(l-y) 



from which we infer that the n th moment of the endpoint region is given in 
terms of the matrix element {B(v)\h v (iD + ) n h v \B(v )). 

The function / is a universal distribution function, which appears in 
all heavy-to-light inclusive decays; another example is the decay B — »■ X s ^ 
[|H] , |12| , where this function determines the photon-energy spectrum in a 
region of order 1 / uiq around the K* peak. 

In principle / has to be determined by other methods than the I/tuq 
expansion, e.g. from lattice calculations or from a model, or it has to be 
determined from experiment by measuring the photon spectrum in B — > 
or the lepton spectrum in B — > X u iv. In the context of the model ACCMM 
model |22| / has been calculated in p4 . 



Some of the properties of / are known. Its support is — oo < k + < A, 
it is normalized to unity, and its first moment vanishes. Its second moment 
is given by a%, and its third moment has been estimated [12 , 18| . A one- 
parameter model for / has been suggested in [[13], which incorporates the 
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1 1.5 2 

Ei [GeV] 

Figure 2: Charged-lepton spectrum in B —>■ X u iv decays. The solid line 
is fl32|) with the ansatz (^), the dashed line shows the prediction of the 
free-quark decay model. The figure is from |I3|. 



known features of / 
f(k+) = 



32 



7T 



x 



exp|-^(l-x) 2 |e(l-x), 



(38) 



where x = k + /A, and the choice A = 570 MeV yields reasonable values for 
the moments. In fig. |2| we show the spectrum for B — > X u £u£ using the ansatz 



Including the non-perturbative effects yields a reasonably behaved spec- 
trum in the endpoint region and the 5-function-like singularities have disap- 
peared. Furthermore, the spectrum now extends beyond the parton model 
endpoint; it is shifted from E™ ax = itlq/^ to the physical endpoint E™ ax = 
Mjj/2, since / is no n- vanishing for positive values of k + < A = Mjj — itlq. 



6 Conclusions and Open Questions 

The 1/rriQ expansion obtained from the OPE and HQET offers the unique 
possibility to calculate the transition rates for inclusive decays in a QCD 
based and model independent framework. The leading term of this expan- 
sion is always the free quark decay, and the first non-trivial corrections are 
in general given by the mean kinetic energy of the heavy quark inside the 
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heavy hadron Ai and the matrix element A2 of the chromomagnetic moment 
operator. 

The method also allows us to calculate differential distributions, such 
as the charged lepton energy spectrum in inclusive semileptonic decays of 
heavy hadrons. For this case, the expansion parameter is the inverse of 
the energy release — 2E#, where E% is the lepton energy. Close to the 
endpoint, the energy release is small and thus the expansion in its inverse 
powers becomes useless. In this kinematic region one may partially resume 
the 1/tuq expansion, obtaining a result closely analogous to the leading twist 
term in deep inelastic scattering. Particularly in the endpoint region a non- 
perturbative function is needed which corresponds to the parton distributions 
parametrizing the deep inelastic scattering. 

The main focus of this review have been the non-perturbative correc- 
tions, and we did not consider the perturbative ones. In general, the non- 
perturbative contributions are typically a few percent, while the perturbative 
ones are usually in excess of ten percent. An exception is the endpoint re- 
gion of semileptonic decay spectra, where both the perturbative and the non- 
perturbative corrections become large and it becomes hard to disentangle the 
two. 

The method of the I/tuq expansion allows us to systematically calculate 
even purely hadronic inclusive rates. This is remarkable, since calculations 
of exclusive hadronic processes are still very model dependent and thus in 
general not reliable. Having control over the purely hadronic inclusive widths 
means that one is in the position to calculate lifetimes and branching fractions 
within the 1/tuq expansion. 

The leading term of the 1/mg expansion is the free quark decay, and it 
is known that the semileptonic branching fraction calculated in the parton 
model is too large by a few percent. The first non-perturbative corrections 
turn out to be much too small to explain the low semileptonic branching 
fraction. The perturbative corrections are larger and indeed tend to lower 
the semileptonic branching fraction, but the effect is still too small to explain 
the data. Recent calculations IE51 EBI indicate that there are some effects 



originating from the finite mass of the final state quarks, in particular in 
the channel b — > ccs. However, enlarging this channel relative to b — > cud 
will increase the average number n c of charm quarks produced per b decay. 
Experimentally, n c is close to unity, n c = 1.04±0.07 ||27|| , while an explanation 
of the low semileptonic branching fraction via an enhancement of the channel 
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b — > ccs would lead to n c ~ 1.3. 

The problem of the semileptonic branching fraction is at the level of a 
two standard deviation discrepancy, and another problem of about the same 
significance are the lifetimes of b hadrons. Lifetime differences between B 
mesons should show up a the level of the 1 / m\ corrections and thus should be 
small. This is supported by data; however, the situation is different for the 
lifetime of b baryons. Their lifetime may differ at the 1 / ml level and hence 
are expected to be of the order of a few percent. Experimentally one finds 
a large lifetime difference between the A& and the B mesons: t(A&)/t(.BA = 
0.72 ±0.11. 

If with new and improved data these two problems persist, they will need 
clarification and perhaps will lead to new insights into the strong interaction 
aspects of heavy flavour weak decays. 
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